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Description 

METHOD AND PRODUCT FOR 
PROCESSING DIGITAL IMAGES 

Background of Invention 

[0001] The present disclosure relates generally to a method and 
product for processing digital images, and particularly to 
a method and product for processing digital images with 
detection and suppression of background noise. 

[0002] a digital image is a collection of digital information in the 
form of pixels that may be processed to provide a visual 
image representative of an object-of-i merest. The digital 
image may be obtained from digital cameras, digital 
video, digital scanners, or the like, or may be digitized 
from an analog film. An exemplary object-of-interest is a 
biological object, and an exemplary digital scanner is a 
magnetic resonance imaging (MRI) scanner. While refer- 
ence may be made herein to MRI as an exemplary scan- 
ning method, it will be appreciated that other scanning 
methods may be employed, such as computed tomogra- 



phy (CT), ultrasound, or X-ray, for example. Digital im- 
ages of objects that are smaller than the field of view of 
the scan contain background regions where the signal in- 
tensity should be zero. However, variations in the intensi- 
ties of these background pixels may arise due to noise 
and other artifacts such as ghosting and gradient warping 
correction for example. Though noise may be prominent 
throughout an image, it is most prominent in the back- 
ground regions where the signal is expected to be zero. 
The presence of background intensity in an image unnec- 
essarily distracts an observer who may be viewing the fi- 
nal image. Accordingly, there is a need in the art for a 
digital image processing method and product that im- 
proves image quality by detecting and suppressing back- 
ground noise. 
Summary of Invention 

[0003] Embodiments of the invention include a method for pro- 
cessing a digital image. A foreground region relating to an 
imaged object is estimated, a background region relating 
to other than the imaged object is estimated, and by using 
the image, the estimated foreground region and the esti- 
mated background region, a transition region disposed 
between the foreground region and the background re- 



gion is calculated. The estimated foreground region, the 
estimated background region, and the calculated transi- 
tion region, each include a separate set of pixels that may 
each be processed separately for suppressing pixel inten- 
sities in the estimated background region and improving 
image quality. 

[0004] other embodiments of the invention include a computer 
program product for processing a digital image, the prod- 
uct including a storage medium, readable by a processing 
circuit, storing instructions for execution by the process- 
ing circuit for performing embodiments of the aforemen- 
tioned method. 
Brief Description of Drawings 

[0005] Referring to the exemplary drawings wherein like ele- 
ments are numbered alike in the accompanying Figures: 

[0006] Figure 1 depicts an exemplary method and system for 
processing digital images in accordance with embodi- 
ments of the invention; 

[0007] Figure 2 depicts further detail of the method of Figure 1; 

[0008] Figure 3 depicts further detail of the method of Figure 2; 

[0009] Figure 4 depicts an alternative method to that depicted in 
Figure 3; and 



[0010] Figure 5 depicts yet further detail of the method of Figure 
1. 

Detailed Description 

[0011] Embodiments of the invention provide a method and 

product for processing digital images, and particularly to 
a method and product for processing digital images with 
detection and suppression of background noise. In an ex- 
emplary embodiment, a biological object is scanned using 
an MRI imaging technique, and the resulting image is dig- 
itally processed for the detection of background noise and 
suppression thereof, thereby improving the signal- 
to-noise-ratio (SNR) of the resulting image. The digital 
image is separated into foreground, background and 
transition regions, thereby enabling each region to be an- 
alyzed and processed separately. The processing tech- 
niques applied involve low-level image processing tech- 
niques, such as thresholding and binary image subtraction 
for example. While embodiments described herein may 
employ MRI as an exemplary imaging technique, it will be 
appreciated that the disclosed invention may also employ 
other imaging techniques, such as CT, Ultrasound, X-ray, 
and the like, for example. 

[0012] Referring now to Figure 1, a high-level flowchart of an ex- 



emplary method 100 for processing a digital image is de- 
picted. An original image 110 is analyzed and segmented 
at block 120 into foreground, background and transition 
regions, with the background region being referred to as a 
background mask 130, and the foreground and transition 
regions being incorporated into and referred to as a tran- 
sition mask 140. A background detection method 200, 
discussed hereinafter in reference to Figure 2, segments 
original image 110 into the three aforementioned regions. 
While reference is made herein to an "original" image 110, 
it will be appreciated that this reference is merely a refer- 
ence to the image entering analysis block 120, and does 
not infer a "first occurrence" of image 110. For example, it 
is contemplated that original image 110 may be a digital 
representation of an analog film resulting from an imag- 
ing technique, such as X-ray for example. Here, the ana- 
log film would be considered the "first occurrence" image, 
and the digital representation would be considered the 
original image 110 entering analysis block 120. 
[0013] information from original image 110, background mask 
130 and transition mask 140, is used by a background 
noise suppression method 300, discussed hereinafter in 
reference to Figure 3, for suppressing the background 



noise and improving the SNR of the resulting image. The 
resulting image is then filtered to produce a filtered image 
150. While embodiments described herein disclose the 
suppression of background noise, it will be appreciated 
that such suppression may generally be described as the 
suppression of pixel intensities in the background region. 
[0014] Background detection method 200 will now be discussed 
with reference to Figures 2-4. From original image 110, a 
foreground region relating to the imaged biological object 
is estimated 202 by applying a pixel intensity threshold 
t_high 402 to filter out all but those pixels having an in- 
tensity equal to or greater than t.high 402, where such 
filtering may alternatively be viewed as filtering original 
image 110 according to a first criterion. In an exemplary 
embodiment, the value of t_high 402 is set to be about 
50% of the mean intensity of the non-zero pixels in origi- 
nal image 110. The resulting region is referred to as the 
estimated foreground region 204, also referred to as the 
initial foreground region. Also from original image 110, a 
background region relating to other than the imaged ob- 
ject is estimated as containing all pixels in original image 
110 having an intensity of zero. The resulting region is 
referred to as the estimated background region 206, 



which is also referred to as the zero mask. Alternatively, 
the estimated background region 206 may be defined as 
that region containing those pixels of the image meeting 
a second criterion, such as an intensity of zero, or an in- 
tensity less than a low intensity threshold, for example. 
[0015] From original image 110, a gradient magnitude image 208 
is computed that provides a value for the gradient of the 
intensity of each pixel of original image 110. A gradient- 
constrained-hysteresis-threshold (GCHT) method 400 is 
applied to original image 110 using gradient magnitude 
image 208 to calculate an initial transition region 210. 
Initial transition region 210 is disposed between esti- 
mated foreground and background regions 204, 206, 
which may be seen by referring to the illustration depicted 
Figure 5. 

[0016] | n an embodiment, initial transition region 210 is calcu- 
lated to be that region containing pixels having a mor- 
phological connection (that is, a connection made possible 
via a morphological operation) to a pixel of estimated 
foreground region 204, having an intensity greater than a 
low threshold t_low 404, and having a gradient magnitude 
that is within a gradient tolerance value g_tol 406 of the 
gradient magnitude of the foreground pixel to which it is 



connected, which is depicted in algorithm form in Figures 
3 and 4. Generally, however, initial transition region 210 
may be considered to be that region of the image con- 
taining those pixels that meet a transition region criterion, 
such as the set of aforementioned criteria, for example. 

[0017] Referring to Figure 3, the application of GCHT method 

400 with hysteresis thresholding and gradient constrain- 
ing is depicted in the absence of an iterative approach. 
Here, the inputs to GCHT method 400 are gradient mag- 
nitude image 208 and estimated foreground region 204, 
where estimated foreground region 204 has been esti- 
mated using intensity threshold t_high 402. The output of 
GCHT method 400 is a gradient constrained hysteresis 
threshold image 450. For each pixel x in the foreground 
(estimated foreground region 204 for example) where the 
pixel intensity Ix is greater than t_high 402, initial transi- 
tion region 210 is calculated to be that region defined by 
the set of pixels y where: 

[0018] y j S connected to x in an 8-connected sense (that is, in a 
3-by-3 matrix of 9 pixels, the center pixel is considered 
to be connected to each of the 8 pixels on the perimeter); 

[0019] the intensity ly of y pixel is greater than low threshold 
t_low 404; and 



[0020] the absolute value of the difference between the gradient 
magnitude of x pixel and the gradient magnitude of y 
pixel (|Gx-Gy|) is less than a defined gradient tolerance 
value g_tol 406. 

[° 021 ] In GCHT 400, the values for low threshold tjow 404 and 
gradient tolerance value g_tol 406 may be user adjusted, 
thereby making GCHT method 400 tunable. 

[0022] Estimated foreground region 204, estimated background 
region 206, and calculated initial transition region 210, 
each comprise a separate set of pixels that may each be 
processed separately for suppressing background noise 
and improving image quality. For example, and as sug- 
gested by GCHT method 400 depicted in Figure 3, the y 
pixels of initial transition region 210 may be added to the 
x pixels of estimated foreground region 204, thereby 
growing the foreground region that is representative of 
the biological object being imaged for greater object ren- 
ditioning. Here, and in the absence of an iterative ap- 
proach, initial transition region 210 becomes transition 
mask 140. 

[0023] An advantage associated with GCHT method 400 is that in 
an image having strong edges, the growing region grows 
along the edge and not across the edge. In this way, mod- 



erate intensity background artifacts that are connected to 
an edge will not get incorporated into the transition re- 
gion. As a result, the t_low threshold 404 may be set quite 
low to capture low intensity foreground regions into the 
transition region without incorporating artifacts with 
higher intensities that are in the background. Also, this 
approach will capture the low intensity foreground regions 
in images that suffer from intensity variation as long as 
the image gradient in the regions of intensity variation is 
less than the gradient tolerance value g_tol. 

[0024] Referring to Figure 4, GCHT method 400 is modified to 

include an iterative approach, thereby resulting in an iter- 
ative GCHT method 500. Here, GCHT method 400 is 
shown as one operation in an iterative do-loop 505. 

[0025] At the beginning of the iterative process, the gradient tol- 
erance value g_tol 406 is set 510 to a low gradient toler- 
ance value toLlow 408, after which, GCHT method 400 is 
applied as discussed previously. Next, at 520, it is deter- 
mined whether the foreground-plus-transition-region (for 
example, estimated foreground region 204 plus initial 
transition region 210 for the first pass in the iterative pro- 
cess) has grown by more than a defined number of pixels 
size_t. In response to the determination at 520 being yes, 



all pixels that may be added to the foreground since the 
last labeling are labeled as incremental transition region 
"i", and the iteration counter is incremented by one 
(i = i + l). These actions are depicted at 530. In response to 
the determination at 520 being no, gradient tolerance 
value g_tol 406 is incremented by a gradient tolerance 
value step g_tol_step 410 until a high gradient tolerance 
value toLhigh 412 condition is met or exceeded. These 
actions are depicted at 540 and 550, respectively, where 
550 is labeled "repeat while: g_tol<tol_high". With each it- 
eration of do-loop 505, incremental transitions regions (i, 
i + 1, i+2, etc.) are calculated, with each incremental tran- 
sition region having an incrementally larger gradient tol- 
erance value g_tol 406 that has a value between toljow 
408 and toLhigh 412. The iterative do-loop 505 is con- 
tinued until toLhigh 412 is met or exceeded. 
[0026] | n an embodiment employing the iterative method of Fig- 
ure 4, each incremental transition region (i, i + 1, i+2, etc.) 
530 is calculated to be that region containing pixels con- 
nected to a pixel of a previously calculated incremental 
transition region 530, having an intensity greater than a 
low threshold t_low 404, and having a gradient magnitude 
that is within the incrementally larger gradient tolerance 



value (g_tol+g_tol_step) of the gradient magnitude of the 
incremental transition region pixel to which it is con- 
nected. 

[0027] upon completion of do-loop 505, that is, where the con- 
dition "repeat while:g_tol<tol_high" 550 is no longer sat- 
isfied, all of the incremental transition regions (i, i + 1, i+2, 
etc.) 530 are merged into a single transition region 560 
using a focus parameter 570. In an embodiment, the fo- 
cus parameter 570 is the percentage of the total number 
of transition regions that will be kept as transition mask 
140. For example, where there are ten transition regions 
(i= 1 through 10) and the focus parameter 570 is set to 
80% (focus=0.8), then the first eight transitions regions 
(i = l through 8) will be merged into transition mask 140 
and the last two transition regions will be discarded. Ac- 
cordingly, and in the presence of an iterative approach, 
the adjustment of single transition region 560 by focus 
parameter 570 results in transition mask 140. An advan- 
tage of using a focus parameter is that the internal pa- 
rameters of method 100 may be fixed and the entire pro- 
cess may be made controllable or tunable by the single 
focus parameter 570. The higher the focus parameter 
570, the more liberal method 100 is in accepting pixels 



into transition mask 140. 
[0028] | n an embodiment, the following parameters may be tun- 
able: 

[0029] intensity high threshold tjiigh 402 

[0030] intensity low threshold tjow 404 

[0031] gradient tolerance value g_tol 406 

[0032] | ow gradient tolerance value toljow 408 

[0033] gradient tolerance value step g_tol_step 410 

[0034] high gradient tolerance value toLhigh 412. 

[0035] However, some of these parameters may set to a defined 
value. For example, low gradient tolerance value toljow 
408 and high gradient tolerance value toLhigh 412 may 
be set to a specific percentage of the average gradient 
magnitude of the entire original image 110, and gradient 
tolerance value step g_tol_step 410 may be set according 
to a desired number of iterations, such as 
(tol_high-tol_low)/10 where ten iterations are desired for 
generating ten incremental transition regions 530, for ex- 
ample. 

[0036] Once transition mask 140 is defined, an object region, 
also referred to as the object image, may be defined as 



being the union of the estimated foreground region 204 
and the transition mask 140, which may then undergo 
morphological operations for improving the image quality. 
Here, the term object region refers to those pixels of the 
original image that are determined via method 100 to be 
representative of the biological object under observation. 
Such morphological operations may include erosion of the 
image and dilation back to the transition mask to remove 
small objects and bridges that may connect background 
noise to the main object, and dilation of the image and 
erosion back to the original size to fill small cracks in the 
edges of the image 212. The zero mask may then be 
combined with the object image, which may then undergo 
a connected-components morphological operation to fill 
small holes in the resulting image 214. Upon completion 
of the hole filling process 214, the zero mask is removed 
from the image resulting in the final object mask 216. The 
final object mask 216 is then separated into a final fore- 
ground region and a final transition region 218. In accor- 
dance with embodiments of the invention disclosed 
herein, the initial estimated foreground region 204 is la- 
beled the final foreground mask 152, the difference be- 
tween the object region and the final foreground mask is 



labeled the final transition mask 154, and the remainder 
of the image is labeled the final background mask 156, 
which collectively make up the filtered image 150, best 
seen by referring back to Figure 1. By separating the im- 
age into different regions, each region may be treated dif- 
ferently to improve image quality. For example, the inten- 
sities of the foreground pixels may be increased slightly 
to improve contrast and increase SNR in the foreground 
regions, while the intensities of the background pixels are 
lowered to reduce noise in the background regions. The 
final object image is the union of the final foreground 
mask 152 and the final transition mask 154 after filtering. 
[0037] Regarding the suppression of background noise, and re- 
ferring now to Figure 5, the background noise may be 
suppressed gradually by reducing the intensity of the 
background pixels to zero as a function of their distance 
from the final object image. In Figure 5, final foreground 
mask 152, final transition mask 154, and final back- 
ground mask 156, are depicted. Shown at increasing dis- 
tances from final transition mask 154 are two regions 
310, 320, which extend into the background region de- 
fined by final background mask 156. In accordance with 
embodiments of the invention, the background noise in 



regions 310, 320 is suppressed gradually by reducing the 
intensity of the background pixels in those regions to 
zero the further the region is from final transition mask 
154. Different functions may be employed to gradually re- 
duce the intensity of the background pixels to zero, such 
as a Gaussian function, a Hanning function, a Hamming 
function, or any function suitable for reducing a pixel in- 
tensity value with respect to distance from a defined ori- 
gin. Combinations of functions may be employed as de- 
sired. Alternatively, all background pixel intensities may 
be identically set to zero, or the background region may 
be homogenized by iteratively blurring the background 
region and adding a small fraction of the original image 
back into the background region. 
[0038] Embodiments of the invention may be provided in exe- 
cutable instruction form on a storage medium, such as 
memory 605 or in the form of a hard drive, a floppy disk 
or a CD-ROM 610 for example, that is readable by a pro- 
cessing circuit, such as processor 615 for example, the 
processing circuit being in signal communication via ap- 
plication software with a graphical user interface at a 
computer, such as computer 600 for example, whereby a 
user may execute the embedded instructions for practic- 



ing the disclosed invention. The instructions may be 
loaded into and/or executed by computer 600, or trans- 
mitted over some transmission medium, such as over 
electrical wiring or cabling, through fiber optics, or via 
electromagnetic radiation, wherein when the computer 
program code is loaded into and executed by computer 
600, the computer becomes an apparatus for practicing 
the invention. When implemented on a general-purpose 
microprocessor, the computer program code segments 
configure the microprocessor to create specific logic cir- 
cuits. The technical effect of the executable instructions is 
to enhance the quality of a digital image, of a scanned bi- 
ological object for example, by detecting and suppressing 
background noise in the image. 
[0039] As disclosed, some embodiments of the invention may in- 
clude some of the following advantages: use of low level 
image processing tools (no statistical models, no neural 
network based methods, for example) for improving im- 
age quality; use of method parameters that may be tun- 
able or predefined; the ability to capture low intensity 
foreground regions into the transition mask without in- 
corporating artifacts with higher intensities that are in the 
background; and, the ability to execute morphological and 



filtering operations on discrete regions of the image 
thereby providing greater control over the generation and 
quality of the final object image. 
[0040] while the invention has been described with reference to 
exemplary embodiments, it will be understood by those 
skilled in the art that various changes may be made and 
equivalents may be substituted for elements thereof with- 
out departing from the scope of the invention. In addition, 
many modifications may be made to adapt a particular 
situation or material to the teachings of the invention 
without departing from the essential scope thereof. 
Therefore, it is intended that the invention not be limited 
to a particular embodiment disclosed as the best or only 
mode contemplated for carrying out this invention, but 
that the invention will include all embodiments falling 
within the scope of the appended claims. Moreover, the 
use of the terms first, second, etc. do not denote any or- 
der or importance, but rather the terms first, second, etc. 
are used to distinguish one element from another. Fur- 
thermore, the use of the terms a, an, etc. do not denote a 
limitation of quantity, but rather denote the presence of at 
least one of the referenced item. 



